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^ Abstract 

We study the ground state entanglement properties of a finite size ensemble of interacting qubits 
driven by a quantized field. In the absence of driving, the system reduces to a Lipkin-Meshkov-Glick 
model that presents maximal entanglement at the first of a series of quantum phase transitions; i.e. 
a 14^-state of the qubit ensemble. We present exact and numerical solutions, within and without 
the rotating wave approximation, in that order, showing that a large finite probability of having 
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the maximal entangled state in the ground state is preserved in the driven system. 
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Entanglement is a fundamental quantum mechanical phenomenon [T] and a precious re- 
source in quantum information processing [2J. A multi-qubit W-state [3] maximizes pairwise 
entanglement of formation [U[5] and is a robust source of entanglement [SIC]; i-e., it retains 
maximal bipartite quantum correlations whenever any pair of qubits are traced out. Interac- 
tion between qubits is necessary to generate entangled states. For an ensemble of interacting 
qubits, the Lipkin-Meshkov-Glick (LMG) model is an example of a system producing max- 
imal pairwise entanglement of formation at the phase transition of its ground state; e.g., 
a collection of spins with (anti-)ferromagnetic interactions embedded in a magnetic field 
produce maximal entanglement at the (first)second order quantum phase transition [8]. 

Recently, it has been proposed that the LMG model, originally from nuclear physics, 
may describe the Josephson effect or two-mode Bose-Einstein condensate (BEC) [9]. Also, 
schemes to realize a dissipative LMG model in optical cavity-QED pU] and in circuit- 
QED [TJ] have been discussed. 

In this work, we consider a system consisting of an ensemble of interacting qubits driven 
by a quantized field and described by the extended Dicke Hamiltonian, 

H = u f N + 6J Z + vK 1 ^ + w~ 1/2 (a + a f ) j x , (i) 

which is equivalent to the LMG model in the limit A = 0. The detuning 5 has been defined 
as the difference between the qubit transition frequency u and the quantized field frequency 
cjf, 5 = oj — cuf, the orbital angular momentum representation, Jj with i = x,y,z,±, has 
been chosen to describe the qubit ensemble, while the field is described by its creation 
(annihilation) operators, a (a)) and the total excitation number is given by = a) a + J z ; 
the field-ensemble and the inter-qubit coupling are represented by the constants A and r\, in 
that order. 

This Dicke-LMG model, Eq.Q, may describe an experimentally feasible two-hyperfine- 
structure-defined-modes BEC coupled to a quantized cavity field mode through a one mi- 
crowave photon process; e.g., trapped hyperfine ground states of a Sodium BEC p2] inside 
a microwave cavity. It may also be possible to realize this model with arrays of interacting 
superconducting qubits [13] coupled to the quantized field mode of a coplanar waveguide 
resonator. 

The ground state phase transition of the Dicke-LMG model has been studied in the ther- 
modynamic limit, iV a — > oo, within and without the rotating wave approximation (RWA) 
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[H] and in the quantum regime, using coherent states for both the field and ensemble, with- 
out the RWA [15] ; these results show the existence of a finite size first order quantum phase 
transition and a second order super-radiant phase transition. Quantum phase transitions 
in the Dicke model have been associated with entanglement between the ensemble and the 
quantized field jT6j H7j and bipartite entanglement among qubits due to finite-size effects 

P2H2U]. 

In the following, we will show analytically, within the RWA, that the maximal entan- 
glement, produced by the quantum phase transition, in the Lipkin-Meshkov-Glick model 
is retained, in the weak field-ensemble coupling regime, by the interacting qubit ensemble 
coupled to the quantized field. Furthermore, for a fixed field-ensemble coupling, A, the inter- 
qubit coupling, rj, defines a series of first order phase transitions related to the number of 
qubits in the ensemble, N a , as the quotient N a /2, while, for a fixed inter-qubit coupling, 
there exists one second order phase transition related to the couplings ratio [15J and a series 
of first order transitions similar to those in the Dicke model jTHJ [19]. Also, we will show 
analitically in the weak coupling regime, and numerically in general, that the inclusion of 
counter-rotating terms does not erase the possibility of obtaining a maximally entangled 
multi-qubit state. 

The ground state of Eq.([T]) has been approximated using the Holstein-Primakoff trans- 
formation and the boson expansion method. It was shown that a normal to superradiant 
phase transition occurs at zero temperature [14J; also, boson coherent states for the field 
and the Dicke basis for the qubit ensemble (T5] have recently been used to show the ex- 
istence of both a first order quantum phase transition and a second order superradiant 
phase transition. However, it is possible to find the exact ground state for the Hamilto- 
nian in Eq.([T]) in the RWA by partitioning the corresponding Hilbert space into subspaces 
with a mean total excitation number (N) = n — N a /2 for n = 0,1,2,... due to con- 
servation of the total number of excitations, [H, N] = 0. The ground state at the n-th 
subspace is given by \ipg^) = Y^k=o c aT' > l^)/l n — k — N a /2) such that H\ipg) = Eg n ^\i/j g n ^) 
and Eg is the lowest eigenvalue for the given subspace. The ket notation |n)/ refers to 
a Fock state of the field with n photons and the shorthand notation |m) = \N a /2,m), 
m = —N a /2, 1 — N a /2, . . . , N a /2 — 1, N a /2, is used for the Dicke state corresponding to the 
normalized superposition of all possible combinations of N a atoms with (N a /2 + m) of them 
in the excited state and the rest, (N a /2 — m), in the ground state. The ground state is found 
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as the lowest eigenvalue for the set {H^ = Hq +Hj n ^} of square matrices, with rank equal 
or less than (N a + 1), defined by the auxiliary matrices, 

=u f {n-N a /2)l fl , (2) 
(H[ n) ) = Sijdt + AiV- 1/2 (Sij-KJj + 5 id+1 o j+1 ) , (3) 

where the identity matrix of rank d is given by Id, the row and column labels are in the range 
i, j — h, h + 1, . . . , n, where n = max(0, n — N a ), for the photon number n — 0, 1,2, . . . The 
symbol stands for Kronecker delta and the diagonal and off-diagonal terms are defined 

as, 

d . = ( n -j- NJ2) [6 + V (n - j - N a /2)/N a ] , (4) 
o, = [j(N a +j-n)(n-j + l)} 1/2 . (5) 

The auxiliary matrix is a tri-diagonal symmetric real matrix with positive off-diagonal 
terms, i.e., a Jacobi matrix, and its eigenvalues can be found analytically [2"Tl - [2~4"] . Further- 
more, at the intersection of two ground state energies belonging to contiguous subspaces, a 
first order quantum phase transition is located. 

The first ground state structure, which will be called vacuum phase from now on, corre- 
sponds to the vacuum field and the atomic ensemble ground state, \ipg°^) = |0)| — N a /2). A 
first quantum phase transition, in a series of first order quantum phase transitions, is found 
at the critical coupling strength, 
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X cl = {[oo + (l/N a -l) V ]u f }2 , (6) 
with the corresponding first non-vacuum state, 

IVf) = 4 1} |0) |1 - NJ2) + c«|l) |-JV a /2> • (7) 

The amplitudes are given by the expressions Cq = h/(h 2 + I) 1 ' 2 and cS 1} = l/(h 2 + I) 1 ' 2 
related to the amplitude parameter, 

h = (1 - l/N a ) V + 5 - {AX 2 + [(1 - l/N a ) V - 5] 2 Y' 2 ^ 

2A 

Any extended Dicke model that conserves the total number of excitations has a first non- 
vacuum phase of the form given by Eq.(j7|), which includes the ground state of the atomic 
ensemble in the maximally entangled W-state |1 — N a /2). The studied model, in the limit 



h — y yields amplitudes Cg — )■ and — )■ 1; z.e., the state is the fully separable state 
|1) \—N a /2). In the opposite case, h — > oo, the amplitudes are — > 1 and — > and 
the state is the product of the vacuum field and the maximally entangled ensemble W-state, 
|0) |1 — N a /2). In the intermediate cases, h G (0, oo), there is a finite probability of finding 
the ground state for this first non-vacuum phase in the maximally entangled W-st&te. Note 
that the case h — > oo can be obtained for small coupling strengths, A — > 0; i.e., the maximally 
entangled ground-state in the studied model is produced by the intra-ensemble interaction 
as expected by the fact that null field-ensemble coupling, A = 0, delivers the LMG model. 

In the absence of intra-ensemble interaction, t] = 0, the Dicke model is obtained. It is well 
known that the shared bipartite concurrence at the first of a sequence of first order quantum 
phase transitions in the Dicke model is not maximal, nor is it at the subsequent first order 
quantum phase transitions [18]. Other extended Dicke models also lack the characteristic 
of a maximally entangled ground state at any first order quantum phase transition as well; 
e.g., that includes a quantized field non-linearity [20J. 

In order to find out if the maximally entangled W^-state shown at the first non-vacuum 
phase survives the inclusion of counter-rotating terms, it is possible to use the unitary 
transformation, 

U = e -^ a+at ) J «, f = 2\N- 1/2 /(cu + cuf). (9) 

In the weak coupling regime, A C w such that £ <C 1, and requiring a weak intra-ensemble 
interaction, r\ <C uj, it is possible to approximate, 

Hcr = U~ HcrU, 

w u f a ] a + ooJ z + i]N- l j 2 z + ~\N- 1/2 (aJ + + a f J_), (10) 

where the auxiliary coupling strength is given by A = 2ufX/(u + ojf). This is the original 
Hamiltonian with a modified coupling strength, i. e. the unitary transformation under weak 
coupling and intra-ensemble interaction acts as a rotating wave approximation, which has a 
first order quantum phase transition at the critical value on the weak coupling regime given 
by the expression, 

'CO + (1/JVa - 1)T]~ 1 



{C R) _ (g + Wf) 
A cl - o 



(11) 



Note that on resonance the expression for the first critical coupling in the weak coupling 



regime, Eq.(ll ), is equal to the critical coupling in the rotating wave approximation, Eq.(p 



The ground state at the first non- vacuum phase is described again by Eq.(u| substituting 
A —7- A. The maximally entangled W-staie survives the inclusion of counter-rotating terms 
for the region of interest; i.e., the weak coupling regime. 

In the strong coupling regime, A 3> u>, it is possible to suppose the state of the system as 



the superposition of coherent states, \a)j \9, 
this representation 
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By minimizing the mean value energy for the three coherent parameters, {a,8,<f)}, one 
can find the critical coupling in the strong limit; but this goes beyond the scope of this 
manuscript. 

It is known that finite size and classical limit analyses, i. e. an infinitely large ensemble in- 
teracting with a coherent field, for Hamiltonians related to Dicke model yield identical critical 
couplings [18-20]. In order to round up the given analysis, we calculate the phase transition 



in the large number of particles limit. By using the Holstein-Primakoff transformation [26J, 

J z = b^b — N a /2, J + = N^ 2 & ^1 — b^b/N^j , J + = Jl, and a Taylor expansion of the 
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square root in terms of (1/N a ), up to first order in (1/N a ), (l — ¥b/N a j w 1 — b*b/2N a , 
up to a constant energy term of N a (r)N a — 2u>)/2, the Hamiltonian in Eq.([T|, within the 
RWA, transforms into the expression, 



H C l ~ Ufa'a + (u — T]N a )Pb + r](tfb) 



+A 



iftt 1 - 



2N a 



(15) 



In a Fock basis, the system presents a sequence of first order quantum phase transitions. In 
order to find the first of these, the problem transforms into diagonalizing the Hamiltonian 
in the one excitation subspace spanned by the states {|0, 1) , |1, 0)}. Thus, the ground state 
energy in the classical regime for the first non- vacuum phase is given by the expression, 
1 



^CLg - 2 



+ ou f + (l-N a )r ] -{4\ 2 + [uj-UJ f + (l/N a -l)r ] Yy 2 • (16) 



The first critical coupling can be found for Eg = and yields a form identical to the exact 
first critical coupling strength found above in Eq.([6]), 



ACL) 



{{u + {l/N a -l)rj\u f Y. (17) 



Equation ([I]), by using the Holstein-Primakoff transformation and the expansion of the 
square root the Hamiltonian, as stated above, transforms into the following Hamiltonian, 



l CLCR 



ufa ] a + (u - r])b r b + riN^itib) 2 + A (a + a f ) (S + 6 f ) 
-A(2iV a )- 1/2 (a + a t )S t (b + b^b. (18) 

In order to find the phase transition, it is possible to use the coherent state basis \a) a \(3) b , 
i.e. a mean-field approximation, such that the mean value for the energy is given by 



(H CR ) « uj f \a\ 2 + (u- V ) \(3\ 2 + r,N- 1 1/3| 4 

+ (A - A(2iV a )- 1 \(3\ 2 ) (a + a*) (/3 + (3*) (19) 



Optimizing with respect to the coherent parameters a and (3, assuming them to be real 
numbers for the sake of simplicity, leads to a trivial solution and a non-trivial solution that 
requires A > [(u — rfujf] 1 / 2 /2, thus a critical coupling strength for the Hamiltonian including 
counter rotating terms in the classical regime is found, as 

This critical coupling, X^ LCR \ at which a second order superradiant phase transition occurs 
[To] , is half of the critical strength found for the case without counter-rotating terms, A^ , 
as (1/JV ) ~ for ^> 2. This is similar to what happens for the Dicke model in the 
classical limit result where accounting for counter-rotating terms halves the critical coupling 
found without the counter-rotating terms [27] . Also, this last result, X^ L , is directly 
related to that in the literature p~U US] . 



In order to verify our analytical results, in Fig.([T]), we show a numerical phase diagram 
defined by the coupling strength, A, and interaction parameter, 7], for a finite size qubit 
ensemble, N a = 5. The maximum shared bipartite concurrence following the entangled web 
approach jl] and the field-ensemble entanglement probed through von Neumann entropy, 
(S) = — Tr a (p a log pa), also known as entropy of entanglement [2], of the reduced qubit 
ensemble obtained after tracing out the field from the total density matrix, p a = tr/p, are 
shown in Fig. ([I]) for the case within the RWA and in Fig. ^ for the full model. The first 
and second critical couplings shown in dashed black lines in Fig.Q are analytically exact 
obtained from the analysis presented above. As analytical and numerically exact results 
agree for this case, only the numerical graphics are shown. Likewise, in Fig.(|2]), by following 
an equivalent approach to that presented in reference [20], numerical results and the first 
critical coupling including counter-rotating terms are shown. Note that, for this case, it is 
not possible to obtain closed form exact results but the analytical critical coupling found for 
the weak coupling regime is in agreement with the results obtained numerically. 

A maximally entangled ground state is found for very small values of the coupling 
strength, such that 77/ A 3> 1. Stability analysis of the Hamiltonian in the classical limit 
show that a pitchfork bifurcation point of the classical dynamics may occur at these cou- 
pling strength ratios [2E], as conjectured for coupled giant spins [29|. 

Figures [T|a) and [l](b) show that it is possible to have a VF-ground-state for a penta- 
partite ensemble as expected in the weak coupling limit. For five qubits, there exists a large 
parameter region at the first non-vacuum phase where maximal value for the maximum 
shared bipartite concurrence, Cw = 2/5, is reached. This value wanes subsequently as 
expected. The entropy of entanglement rises to a maximum after the first zone, which is 
related to a second order phase transition, a behavior shared by all studied cases, N a e [2, 6], 
to decay for very large values of the coupling strength, A lo. 

In contrast to the Dicke model, where the first critical coupling is independent of the 
ensemble size, in this extended Dicke model the ensemble size plays a role in the first critical 
coupling, e.g., Eq.([6]). As the involved term is of the form — iV -1 , it stops playing an 
important role for large ensemble sizes. In the Dicke model, the vacuum phase remains fixed 
as the first critical coupling stays constant for an increasing number of atoms, while the 
phases beyond the vacuum phase reduce to allow just the vacuum and classical superradiant 
phases for an infinitely large atomic ensemble. In this extended Dicke model, the area of all 
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the different phases reduces as the ensemble size increases. 

In summary, we have analytically derived the exact critical coupling value and ground 
states at the first quantum phase transition for a BEC-cavity-QED system in the rotat- 
ing wave approximation, including counter-rotating terms, and in the classical limit both 
excluding and including counter-rotating terms. Subsequent ground states of first order 
quantum phase transitions have been shown through analytical derivations and numerical 
simulations. The results obtained correlate among themselves in a similar way as the Dicke 
model. 

While the exact critical coupling at the first of the quantum phase transitions connects 
with known mean-field results, our treatment allows to follow the quantum correlations 
between the field and ensemble as well as among ensemble components. Thus, a striking 
feature is revealed: the existence of a maximally entangled atomic ground state within the 
limits of a certain phase space region at the first non-vacuum phase. The parameter region 
where this maximally entangled W-staie is found decreases as the ensemble grows in size. 
The existence of this maximal entangled qubit ensemble state survives even the addition of 
counter rotating terms. Numerical results for the maximum shared bipartite concurrence, 
following the entangled web approach, and the field-ensemble entanglement, probed through 
von Neuman entropy of the reduced two-level ensemble, were presented. Agreement between 
the exact analytical and the numerical results was obtained. 

We believe that our results show that highly tunable experimental systems realizing 
an interacting ensemble of qubits driven by a quantized field in the weak coupling regime 
may be used as a source of entanglement for quantum information processing; e. g., an ar- 
ray of interacting superconducting flux qubits coupled to a coplanar waveguide resonator 
following the building blocks presented in references [TBI 130] or a two ground-state-hyperfine- 
structure-defined mode BEC driven by a single mode of microwave cavity building on known 
technology like that in Ref . [12] , with the inconvenience that the center of mass motion may 
couple to the internal degrees of freedom [31] . 
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Figure 1. (Color online) The phase diagram of the finite size extended Dicke Hamiltonian for the 
case of a penta-partite qubit ensemble, N a = 5, in the parameter space defined by the field-ensemble 
coupling strength, A, and the intra-ensemble interaction strength, r/. (a) Maximum shared bipartite 
concurrence in the sense of entangled webs, C w . (b) Entropy of entanglement, {S), between the 
field and the qubit ensemble. The corresponding minima and maxima values for the color legend 
are shown below it. The dashed lines show the two first exact critical couplings. 
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